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Abstract 



The Wigner time delay of a classically chaotic quantum system can be expressed semiclassically 
either in terms of pairs of scattering trajectories that enter and leave the system or in terms of 
the periodic orbits trapped inside the system. We show how these two pictures are related on the 
semiclassical level. We start from the semiclassical formula with the scattering trajectories and derive 
from it all terms in the periodic orbit formula for the time delay. The main ingredient in this calculation 
is a new type of correlation between scattering trajectories which is due to trajectories that approach 
the trapped periodic orbits closely. The equivalence between the two pictures is also demonstrated by 
considering correlation functions of the time delay. A corresponding calculation for the conductance 
gives no periodic orbit contributions in leading order. 
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1 Introduction 



Quantum systems whose classical counterparts are chaotic show universal statistical fluctuations that 
are well modelled by Random Matrix Theory (RMT) [1, 2]. In the semiclassical limit, formally as 
h — ► 0, the quantum statistics can be approximated by quantities that involve trajectories of the 
classical dynamics. One focus of work in semiclassics has been to recreate results from RMT by 
semiclassical methods. 

For closed systems, statistical fluctuations in the energy spectrum can be seen, for example, in 
the spectral form factor K{t) which is semiclassically approximated by a double sum over periodic 
orbits. The universal fluctuations from RMT are then due to correlated pairs of periodic orbits. The 
diagonal approximation, comparing an orbit to itself (or its time reverse) leads to the first term in 
the expansion of the form factor [3, 4]. All the other terms in the small time (r < 1) expansion of 
the form factor can be obtained, in agreement with RMT, from periodic orbits with self-encounters. 
These are events where an orbit approaches itself (or its time reverse) very closely so that its partner 
can cross the encounter region differently [5, 6]. Currently it is not known what other type of periodic 
orbit correlations contribute in the regime r > 1, though their contribution is calculated indirectly 
in [7, 8]. Similar methods as for the form factor have been applied, for example, to the conductance 
through open systems where the semiclassical sum is over pairs of open trajectories that start and end 
in the leads [9, 10, 11]. 

For the Wigner time delay in open systems we have the interesting situation that the semiclassical 
approximation can be expressed in two ways. One gives the time delay as a double sum over scattering 
trajectories that enter and leave the system, in a similar way as for the conductance. The other is 
through a relation to a density of states and leads to a semiclassical formula that contains the average 
time delay plus a single sum over the periodic orbits that are trapped in the system. Our main 
motivation for this article is to understand the duality of these two semiclassical pictures. For open 
chaotic cavities, we will start from one of these pictures, the double sum over scattering trajectories, 
and derive all terms in the other semiclassical formula for the time delay. The periodic orbit terms 
are obtained by considering a new type of open trajectory correlation which is linked to the motion 
around periodic orbits. For systems without time-reversal symmetry these periodic orbit encounters 
are sufficient for obtaining the correct periodic orbit terms. For systems with time-reversal symmetry 
one also has to include combinations of self-encounters and periodic orbit encounters. 

We will also consider a correlation function of the time delay. When expressed as a double sum 
over periodic orbits, the diagonal approximation [12, 13] and higher order terms [14] were shown to 
agree with RMT (for a small time expansion) . Using the other semiclassical approximation which is in 
terms of a quadruple sum over open trajectories, it was shown in [15] that the diagonal approximation 
does not give the leading RMT result. This was attributed to the non-unitarity of the semiclassical 
scattering matrix [16]. We will show that the inclusion of off-diagonal terms due to trajectories with 
self-encounters removes the discrepancy with RMT and restores the semiclassical unitarity of the 
scattering matrix. These calculations are similar to those for correlation functions of the conductance 
[17, 11]. 

Our paper is divided as follows. In section 2 we introduce the two semiclassical approaches for 
the time delay. In section 3 we show which correlated pairs of open trajectories recreate the average 
time delay. We then go beyond the average to recreate the periodic orbit contributions in section 4. 
This is achieved by introducing the new type of trajectory correlations which are due to periodic orbit 
encounters. The correlation function of the time delay is considered in section 5, and our conclusions 
follow in section 6. 
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2 The time delay 



For a chaotic cavity with one or more open leads that carry M scattering channels, the incoming and 
outgoing waves are related by the M x M scattering matrix S(E). The Wigner time delay, which 
represents the extra time spent in the scattering process compared to free motion, is defined as [18, 19] 



M-E) = - s Tr 



l^Alndet^^) . (l) 
MdE y ; v ' 



The time delay can be expressed semiclassically both in terms of the trapped set of periodic orbits 
of the open system, and in terms of the open scattering trajectories that enter and exit through the 
leads. 

The description of the Wigner time delay in terms of trapped periodic orbits comes from its relation 
to a density of states which, in general, is the difference between the level density of the open scattering 
system and a free system [20] 

rw(E) = *gw « 2 -£d(E) + , (2) 

Here the density of states d{E) is separated into a mean part d(E) and a fluctuating part d a (E) which 
each have a semiclassical approximation. The approximation for the mean density of states, for a 
chaotic cavity with 2 degrees of freedom, comes from Weyl's law for the corresponding closed system 
d(E) ~ ^l/(2iih) 2 where Q is the phase space volume of the shell of constant energy E. The fluctuating 
part can be expressed, like in the Gutzwiller trace formula [21] as a sum over the periodic orbits. The 
difference is that the sum only includes periodic orbits that are trapped in the system [22, 13]. Using 
these approximations, we can write the time delay as 

Tw( ^ * 1 + i Re £ A pA E)eirSAE) e -^ P (3) 

p,r 

where Th is the Heisenberg time which is related to the average level density by Th = 2irhd(E). The 
first term is the average time spent in the cavity fyy = Th/M. It is equal to the inverse of the classical 
escape rate which can be expressed in the form fi = M/Tn [15]. In the sum p labels the trapped 
primitive periodic orbits and r their repetitions. The orbits have action S p and Maslov index fi p . 
Their stability amplitude A P)T can be written in terms of the stability matrix M p and the period T p 



A p , r = - Tp . (4) 



det(M£ - 
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The description of the time delay in terms of open trajectories comes from the semiclassical approxi- 
mation to the scattering matrix elements [23, 24, 9] 

S ba (E) «-J= £ A^e-'f "« . (5) 

Here S a is the action of the trajectory a and v a is the number of conjugate points along the trajectory 
(plus twice the number of reflections on walls with Dirichlet boundary conditions). The stability 
amplitude A a can be found in [24]. The sum is then over all classical trajectories that start in channel 
a and end in channel b, where the channels fix the absolute value of the angles at which the trajectories 
enter and leave the cavity. From this semiclassical approximation for the scattering matrix elements 
we can obtain an expression for the time delay by substituting into equation (1). When we differentiate 
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the scattering matrix elements we ignore the change in the slowly varying prefactor and only keep the 
term from the oscillating action exponentials 

^-^E E Ta^e^-We-* <"«-V> , (6) 

a,b ce,a'(a— >6) 

where T a = dS a /dE is the time the trajectory a spends inside the system. Here we can see that the 
time delay is a sum over trajectory pairs a, a' both of which start and end in the same channels (a 
and b respectively), followed by a sum over all the possible channels. 

We can also consider a correlation function of scattering matrix elements 

(2,0 

where it is convenient to specify the energy difference in units of Mi^lird)" 1 = hfi, because this will 
simplify the formulae in the following. If we set e = this becomes 

C(0) = Tr \S(E)S^(E)} . (8) 

By using the semiclassical approximation of the matrix elements from equation (5) and expanding 
the action up to first order in energy, S a (E + rj) srf S a (E) + r]T a (E), the correlation function can be 
expressed in terms of pairs of scattering trajectories 

a,fo a,ce'(a—*b) 



from which we can obtain a symmetrized version of the time delay [15] 



= — — Tc 

e =o 2M 



'st (B) MS_S(B)^ 



d£ v ' dE 



(10) 



Equation (10) agrees with the definition of the time delay in (1) because of the unitarity of the 
scattering matrix. If we insert (9) into (10) we obtain again a semiclassical formula for the time delay 
which differs slightly from (6) in that the time T a is replaced by the average time [T a + T a /)/2. Both 
formulae are equivalent and the difference only plays a role in section 5 where it will be discussed. We 
use the relation of the time delay to the function C(e) in the following to simplify the calculation by 
exploiting its similarity to the average conductance of a chaotic ballistic device [10, 11]. In particular 
we will obtain, as in the case of the conductance, simple diagrammatic rules for the semiclassical 
contributions of correlated trajectories. 

Because we can express the time delay equally in terms of open trajectories and trapped periodic 
orbits, there should be a semiclassical equivalence between the two pictures and the following should 
hold 

^E E r tt ^^ei^-^) e -*^-^)«^ + ARe^^ r (^e^^e-^. (11) 

a,b a,a'(a->b) * P,r 

The left hand side is the sum over scattering trajectories, while the right includes an average part and 
a sum over trapped periodic orbits. We shall now show how, by considering the contributions in the 
semiclassical limit of pairs of correlated trajectories from the sum on the left we can recreate all the 
terms on the right. We will consider only systems with two degrees of freedom, but the calculations 
are very similar in higher dimensions [6]. The assumptions for the semiclassical calculations in this 
article are the same as for the conductance in [11]. In particular, corrections related to the finite value 
of the quotient of the Ehrenfest time and the dwell time are neglected. 
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Figure 1: An example of a trajectory (full line) with two encounter regions and its partner trajectory 
(dashed line). The encounter regions are indicated by rectangular boxes. 



3 The average time delay 

In this section we shall derive the average time delay from correlated pairs of open trajectories. For 
systems without time-reversal symmetry the diagonal approximation (pairing a trajectory with itself) 
suffices [15]. For systems with time-reversal symmetry, however, a small correction is needed which 
comes from trajectories that have close self-encounters. The calculation follows similar steps and ideas 
to the calculation of the average conductance [9, 10, 11], which in turn builds on work on spectral 
statistics [6] . We exploit this similarity by concentrating on the correlation function of the scattering 
matrix elements C(e). Its semiclassical approximation is given by (9) which, aside from the exponential 
factor containing the trajectory time T a and a different channel sum, is the same as the conductance 
from [11], so we only highlight the relevant points in the calculation. We shall be more detailed in 
section 4 where we introduce the new correlations that produce the periodic orbit terms. 

The diagonal term for the sum over trajectories that connect channels a and b considers pairs 
where the two trajectories a and a' are identical, and it gives a contribution of 



- T 

a(a— >?>) 



\A a \ 2 e itf * Ta . (12) 



The sum in (12) can be performed by using a sum rule for open trajectories [9] which turns it into an 
integral over the trajectory time T 



oc 



V \A a \\..~ / dTe-^... , (13) 

a(a—>b) 

The exponential term in (13) represents the average probability that a trajectory remains in the system 
for the time T. The sum over channels depends on the symmetry of the dynamics. For systems without 
time-reversal symmetry (re = 1) we can pick both a and b from the M possible channels giving a factor 
of M 2 . For systems with time-reversal symmetry (k = 2) we can also pair the trajectory a with its 
time reversal if the start and end channel are the same (a = b), and this contributes an additional M 
to the channel sum. The diagonal approximation thus becomes 

di M(M + k-1) r {1 _ ie)T = M(M + k-1) 

W T H Jo Af (1 - ie) K ' 

Using (10) this already leads to the correct result for the average time delay for systems without 
time reversal symmetry, but not for systems with time-reversal symmetry (because the prefactor then 
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contains M(M + 1) instead of M 2 ). Hence we need to consider off-diagonal contributions to the 
average of C(e). They come from long trajectories that have one or more self-encounters in which 
two or more stretches of a trajectory are almost identical. In systems with time-reversal symmetry 
the stretches can also be almost mutually time-reversed. The encounter regions are connected to each 
other and to the entry and exit channels by long parts of the trajectory called links. An example of a 
trajectory with two encounter regions is shown in figure 1. The partner trajectory connects the links 
in a different way in the encounter regions, but follows the original trajectory very closely along the 
links. Both encounter regions in the figure are examples of 2-encounters which are traversed by two 
stretches of an orbit. In general an arbitrary number of I > 2 stretches of an orbit can be almost 
identical (up to time-reversal) in an encounter region and one then speaks of an /-encounter. The total 
numbers of the different encounter regions of a trajectory are collected in a vector v whose components 
vi specify the numbers of Z-encounters of a trajectory. For the calculation of the off-diagonal terms 
one has to consider all possible structures or families 3 , i.e. all topologically distinct ways in which two 
trajectories can be correlated. A more precise definition of families can be found in [11], but we shall 
need in the following only the fact that there is a sum rule for the number N(v) of different families 
with the same vector v. Further relevant quantities are the total number of encounters of a trajectory 
V = Yli v h an d the total number of orbit stretches of a trajectory in the encounter regions L = J^i I V[. 
The total number of links of a trajectory is then L + 1. 

The action difference of the trajectories is given in terms of coordinates along the stable and 
unstable manifolds in Poincare surfaces of sections in the V encounter regions. These coordinates 
describe the relative positions of the different stretches of a trajectory. The encounter regions are 
labeled by a and the number of trajectory stretches in them by l a . In the linearized approximation 
the action difference is given by 

V 

S a ~ S a > w ^2 s a u a = su , (15) 

cr=l 

where s a and u a are vectors with dimension (l a — 1), and s and u contain the components of all these 
vectors. 

The summation over the correlated trajectory pairs is simplified by the fact that different families 
with the same vector v give the same contribution [11]. Hence it is convenient to collect all these 
contributions and sum over all trajectories pairs whose correlations are specified by the vector v. As 
for the diagonal approximation, we have also to take into account that one can pair a trajectory with 
the time-reverse of its partner orbit if a = b in systems with time-reversal symmetry, and this gives 
an additional factor of two in these cases. This factor is denoted by 7] a b in the following 

^ fixed v 

C v (e)^jrJ2^ E \A a \ 2 e^- s ^e^ , (16) 

a,b a,a'(a^b) 

where the approximations T a ps T a i, A a ps A a / and fi a pa fi a / have been made, and 

ri ab = 1 + (k - l)5 ab . (17) 

The sum over the trajectory pairs in (16) is performed by applying an ergodicity argument together 
with the finite escape probability of the scattering trajectories which results in replacing it by an 
integral 

fixed v r r 

\A a \ 2 ... rs N(v) dT dsdw^ iT (s,w)e^ Tcxp ... . (18) 

a,a' (a— ►&) 

3 The authors of [11] use the expression structures in the context of periodic orbits and families in the context of open 
trajectories 
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Here N(v) is the number of families with the same vector v, and T is the time of the trajectories. 
T exp will be specified below, and w Vt T(s,u) is the probability density that a trajectory of time T has 
self-encounters specified by the separation coordinates s and u. It is given explicitly in terms of an 
integral over L of the L + 1 link times 

f 1 

w v t(s,u) = / d£i...d£L tj , (19) 

where tj is the time of link i and ig nc is the time of encounter in the encounter region a. The prime 
at the integral denotes that it is subject to the restrictions that the link times must be positive. 
Furthermore, the total time of the links and encounter stretches is the time of the trajectory 



L+l 

£ 

i=i 



ti + X>C:, (20) 



and the time of the last link tL+i, which is fixed by T and the L other link times, has to be positive 
too. The encounter times ig nc are specified by requiring that all the components of s a and u a that 
determine the separation of the different stretches of a trajectory in the encounter region a have a 
modulus that is smaller than a small arbitrary constant c, and it is given by 

Cc ~ t ln r-r^ 1 — i > ( 21 ) 

A maxj \s ai \ x maxj \u aj \ 

where A is the Lyapunov exponent. The relevant encounter times are of the order of the Ehrenfest 
time. Finally, the exposure time T exp differs slightly from the time T of the trajectory, because the 
trajectory stretches during each encounter are very close together. If the trajectory survives during 
one crossing of the encounter region it will survive all the others crossings. The exposure time, the 
effective time where the trajectory can leave, is thus given by T exp = T — ^2 a (la — l)t° nc . Putting 
everything together, i.e. inserting (18), (19) and (20) into (16), one obtains after a change of the 
integration variable T to the last link time tL+i an expression that contains integrals over all link 
times 

CW « -iJ E ,u II (I d ti e-M-.)..) J] [ j d S „ d„ ff ,_ ltgm j . (22) 

One can see that (22) factors into a product over the L + l links and the V encounter regions. This 
is the main advantage of working with the correlation function C(e) instead of directly with the time 
delay, because the corresponding expression for the time delay does not factorize (because of the T a 
in the pre-exponential factor in (6)). The factorization property will be useful also in the following 
sections. The integrals over coordinates in the encounter regions can be performed by using the 
semiclassical result [6, 25] 

d«.d«.(o'.*-~ K {i( 2rt) e.-» HZ?,' < 23 » 

where a is a real constant. The asymptotics of these integrals comes from the origin, or zero separation, 
at which point the semiclassical approximations made above become accurate. We can now perform 
all the integrals in equation (22). The integral over the encounter regions is evaluated by expanding 
the exponential that contains the encounter times up to the linear term. Higher order terms in the 
expansion are neglected, because they are of the order of the quotient of the Ehrenfest time and the 
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dwell time. For every encounter we obtain thus a contribution of — /u(l — iel a )T^~ where /i = M/Th, 
and for every link a contribution of /tx _1 (l — ie) -1 . All the Heisenberg times cancel, and we obtain the 
simple diagrammatic rules that every link contributes a factor of M _1 (l — ie)^ 1 , and every encounter 
a factor of — M(l — id a ). Altogether the result is 

0,6 " v ; 

As for the diagonal approximation the sum over the channels a and 6 gives a factor M(M + k — 1). 
It is convenient to introduce k = L — V + 1 and to combine the contributions from all vectors v with 
the same value of k. We can also include the diagonal term by saying it corresponds to a trajectory 
with no encounter L = V = and N(v) = 1. With the diagonal term, all these trajectories give the 
average contribution to the correlation function 

OO L~V+l = k yrV (-!_■] \ 

(7(e). M(M + .-!)£ £ (-l) y JV(i> fe C _ ffi . (25) 

fc=l w V ' 

We expand the result around e = and use J2a la = L to obtain 

00 L ~V+} =k (_-\\V t\t( \ 

C(e)*M(M + K-l)J2 E [ 1 + ig + °( g2 )]- ( 26 ) 

k=l v 

The sum over the vectors v can now be evaluated by using a sum rule for the numbers N(v) that was 
obtained from a recursion relation in [11] 

L-V+i=k 

E (-l) V iV(«) = (1 - k)*- 1 . (27) 

This leaves the evaluation of a geometric series to obtain 

C{e) m M[l +ie + 0(e 2 )] , (28) 

which is the final result. Equation (28) shows first that the semiclassical approximation is consistent 
with the unitarity of the scattering matrix, because 6(0) = TrS(E)S* (E) ~ M. Secondly, one obtains 
the correct average time delay from the term that is linear in e, fw = Th/M. 

4 Periodic orbit terms 
4.1 Periodic orbit encounters 

The types of correlated trajectory pairs we will consider now approach a trapped unstable periodic 
orbit, follow it very closely, and leave it again. An example of such a trajectory is shown in figure 2. 
The trajectory approaches the periodic orbit almost along its stable manifold, goes around a number 
of times, and then leaves it closely following the unstable manifold. In the vicinity of the periodic orbit 
we describe the motion of the trajectory in a Poincare surface of section transverse to the periodic 
orbit. Because of the Birkhoff-Moser theorem we can make a symplectic transformation to normal 
form coordinates which are along the stable and unstable manifolds. Close to a periodic orbit one can 
use the linearized approximation in which the Poincare map has the simple form [26] 

s^A^s, u' = A p u, (29) 
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Figure 2: A schematic picture of a trajectory that approaches an unstable periodic orbit follows it a 
number of times and leaves it again. In a Poincare map transverse to the periodic orbit the trajectory 
moves along the invariant hyperbola su = const. 

so that the trajectory moves along the invariant hyperbola su = const. A p and 1/A p are the eigenvalues 
of the stability matrix of the primitive periodic orbit p, with |A p | > 1. If A p is negative the map involves 
also a reflection about the origin. The correlated trajectories that we will consider differ in the number 
of times they wind around the periodic orbit, and we shall see that the semiclassical periodic orbit 
contribution is obtained in the limit that su of these orbits goes to zero. This allows the application 
of the linearized approximation. 

We have to specify how to count the number of times a trajectory winds around a periodic orbit. 
This can be done by fixing an arbitrary small positive constant c. The encounter region of the 
trajectory with the periodic orbit is then defined by requiring that the moduli of the coordinates s 
and u are both smaller than c. This is very similar to the definition of the self-encounter regions in 
the previous section. A very long trajectory has a finite probability to enter such an encounter region. 

An example of a trajectory a which has k = 5 intersections with the Poincare surface in the 
encounter region, Pi, . . . ,P$, is shown in figure 3 (for positive A p ). Given such a trajectory one 
can find a partner trajectory with r more intersections in the encounter region in the following way. 
Consider the line through the first point Pi that is formed by the intersections of trajectories that 
satisfy the required initial conditions for scattering trajectories (indicated by the thin line through Pi 
in figure 3). There is second line through the last point Pk due to trajectories that satisfy the required 
final conditions (the thin line through P5 in the figure). If one moves along the first line towards 
the stable manifold then one finds for every r > a unique point P[ such that its (k + r — l)-th 
iterate by the Poincare map, Pfc +r , lies on the second line. This is the required partner trajectory a'. 
Figure 3 shows an example for r = 2. In the following we will show that one obtains the semiclassical 
periodic orbit contributions from this kind of trajectory pairs. We first calculate the action difference 
of the trajectories. As a consequence of the Poincare-Cartan theorem the action integral j pdq is 
independent of the path when it is evaluated on a surface that is formed by a continuous family of 
trajectories on the energy shell [26]. We apply this theorem to calculate three contributions to the 
action difference AS = S a — S a r. AiS coming from the initial parts of the trajectories up to Pi and 
P{, AfS coming from the final parts of the trajectories from P& and P' k+r onwards, and A m S from 
the remaining middle parts 

fPi fP'h+r fO [Pk fO rP' k+r 

AiS = / pdq, A f S = pdq, A m S = / pdq+ / pdq- / pdq - / pdq-rS p . (30) 

JP[ JP k JPx JO JP[ JO 

The first integral is evaluated along the line formed by trajectories with the required initial conditions 
(they are all exponentially close to each other at their starting points) . The second integral is evaluated 
along the line formed by trajectories with the required final conditions (they are all exponentially close 
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Figure 3: Two trajectories which follow the periodic orbit at O five and seven times, respectively, within 
the region of the Poincare section bounded by the constant c. The action difference is S a —S a / = A—2S P 
where S p is the action of the periodic orbit and A is the sum of the areas of the two triangles OP[P\ 
and OP^Plj. 

to each other at their end point). The remaining part of the action difference is obtained by comparing 
the actions of the middle parts of the trajectories to the action of the periodic orbit. The integrals 
are evaluated along lines that are formed by the intersection of trajectories that continuously connect 
the middle parts of the trajectories to, respectively, the /c-fold and {k + r)-fold traversals of the 
primitive periodic orbit p. If we add up all contributions we obtain an integral over the two triangular 
regions OP[P\ and OPkPk+r depicted in figure 3. The areas of these regions are invariant under the 
transformation to normal form coordinates. In the linearized approximation the sides of the triangles 
are straight lines. If P\ = (u, s) then P[ m (uA~ r , s), P k « (uA k , sA~ k ), P' k+r w (uA k , sAp ^ + ^) 5 and 
one obtains 

AS = su(l - A~ r ) - rS p . (31) 

We did the derivation for positive A p , but the coordinates of the points and (31) hold also for negative 
A p . In the case that A p is negative and r is odd then Pi and P[ lie on different sides of the s-axis, 
and Pk and P' k+r on different sides of the n-axis. 

Next we compare the semiclassical amplitudes of the trajectories a and a'. They depend on 
the stability matrices of the trajectories. The stability matrix describes the motion of neighbouring 
trajectories in the linearized approximation. It gives the deviations at the final point of a trajectory 
in terms of the deviations at the initial point, in coordinates perpendicular to the trajectory 

For cavities with leads the amplitude of a scattering trajectory is proportional to l-ZV^il" 1 / 2 where 
M21 is a matrix element of M [9]. The stability matrix can be composed by multiplying stability 
matrices for different parts of the trajectory. For the k iterations of the Poincare map in the vicinity 
of the periodic orbit the stability matrix of a is approximated by the stability matrix of the periodic 
orbit p and hence we write M a = MfM^Mi which i and / stand for the initial and final parts of the 
trajectory. The corresponding approximation for the trajectory a' is M a / = M f M k+r Mi. Powers of 
the two-dimensional stability matrix M p can be written in terms of their eigenvalues as 

M k v = A k p P u + A~ k P s ~ A k p P u as k 00 , (33) 



10 



where P s and P u are the projection operators onto the eigenvectors of M p in the stable and unstable 
directions, respectively. As a consequence M a i ~ M a A p as k — > oo, and we obtain the approximation 

A a , « A a |Ap|- r / 2 . (34) 

It is convenient to have the following geometrical picture. One considers a neighbouring trajectory 
with initial infinitesimal deviation 5qi 7^ and 5pi = and follows the development of this deviation 
in time as one moves along the trajectory. A conjugate point occurs every time the projection of the 
deviation onto the (orthogonal) p-direction is zero, and the index v increases by one. (For billiard with 
Dirichlet boundary conditions it increases also by two for every reflection on the wall.) The matrix 
element M21 is equal to the projection at the final point. During the time the trajectory follows the 
periodic orbit, the deviation is aligned and stretched along the unstable direction, according to (33), 
and rotates with it around the periodic orbit. The Maslov index of the periodic orbit is the number of 
times the stable and unstable manifolds rotate by half a turn [27] (plus twice the number of reflections 
on walls with Dirichlet b.c.s). After each traversal of the periodic orbit, the manifolds are back where 
they started and the Maslov index is an integer. As long as both a and a' are close to the periodic 
orbit, a' will pick up r times the Maslov index of the orbit p over the trajectory a. Outside of the 
encounter, both trajectories are close to each other and have the same index, giving 

Vq! = v a + rn p . (35) 

We mention that this geometrical picture can be generalized to systems with higher degrees of freedom, 
/ > 2. Then one considers the development of an (/ — l)-dimensional volume element in time that 
is formed at the initial point by infinitesimal deviations in the (/ — 1) g-directions orthogonal to the 
trajectory. 

We can now put the trajectory pairs that follow a trapped periodic orbit p in this way into the 
semiclassical sum for C(e), equation (9). The times of the two trajectories are related by T a i « T a +rT p . 
We insert furthermore (31), (34), and (35) into (9) and take into account a factor of two if the channels 
a and b are the same in systems with time reversal symmetry, because then one can pair a trajectory 
also with the time-reverse of its partner. This factor is denoted as before by 77,3,5, see (17). We are left 
with the following to evaluate 

C p ' r (e)^J2^ e " I .a r 2 P £ \A a fe^- A ^e^ + (r^-r). (36) 

The last term denotes the same contribution with r replaced by — r. It is obtained from interchanging 
a and a' in the double sum over trajectories (9). The remaining step consists in performing the sum 
over all trajectories a that enter the encounter region of the periodic orbit p. This summation is 
performed by applying the ergodicity property that very long trajectories explore the available phase 
space volume uniformly, combined with the finite escape probability of the scattering trajectories. As 
a consequence the sum over trajectories, with the amplitudes as weight factors, are replaced by an 
integral (for related sum rules in closed systems see [28]) 

\A a \ 2 ...^ j dT y'dsduu;p i T(s,u)e~' iTcxp ... , (37) 

where T is the total time of the trajectory (=T a ) and w P: t is the probability density that a trajectory 
of time T goes through a surface element ds du around the point (u, s) in a fixed Poincare section 
transverse to the periodic orbit p. The integrals over s and u are limited by the constant c. The time 
T exp is the time during which the orbit can escape. It differs from the time T of the trajectory, because 
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the trajectory cannot escape during the time it follows the trapped periodic orbit in the encounter 
region of the orbit (if one chooses the arbitrary constant c sufficiently small). Hence, the difference 
between T and T exp is the time in the encounter region which is given by 

1 c 2 

— In — , (38) 

Ap | US | 

where k is the number of iterations in the encounter region, and the Lyapunov exponent \ p of the 
periodic orbit is obtained from |A P | = e XpTp . The number k was obtained by using (29) from which 
follows that k « In A + jAjr In A . It is remarkable that (38) has the same form as the encounter 
time for self-encounters (21). For this reason, the central relation (23) for the integrals over s and u 
holds as well and will be used in the following. 

The probability density w Pj t(s,u) is obtained by noting that the probability that an element of a 
trajectory of time dt enters the surface element dsdu is given by 

dsdndt 

n ■ [ ' 

The probability density follows as 

/•T 1 — *enc J 1 

w PtT (s,u)= dti— f-, (40) 

JO * ' Cenc 

where t\ is the time at the intersection with the Poincare surface. Equation (40) contains the factor 
T p /tc nc ~ 1/k to remove the multiple counting of a trajectory, because k points in the Poincare 
section correspond to the same trajectory. We substitute (37) and (40) into (36) and change the 
integration variable from T to t2, where £2 is the time from the encounter region to the exit channel, 
T = t\ + tenc + h, and obtain 

T p-irSp+frfip i^rTp i- jrsu(l-\p r ) p ie/ut§ nc 2 foo 

C^(e) « Y,** — r 71 , £ / ^du G - " J] / d ^ e " Ml " le)tl + (»• - -')■ 

^ T H |A P | 2 y "4 f = Vo 

(41) 

We see again that the expression factorizes into contributions from the links and the encounter region. 
This implies that the diagrammatic rules can be extended, and that one gets additional contributions 
from the periodic orbit encounters. We can now perform the integral over s and u. Note that this 
integral automatically sums over trajectories a with an arbitrary number k of windings around the 
periodic orbit. We expand the integrand in e up to first order and use equation (23). Only the term 
with icnc in the numerator cancels with the encounter time in the denominator and contributes semi- 
classically. After integrating there is again a cancellation of Heisenberg times and the diagrammatic 
rule for the r-th repetition of a periodic orbit p follows as 

2ieM P , r cos (- ~rS p + ^ru p + ^ rT,}j (42) 

where we used the identity yj\ det(M£ — 1)| = |A p | r / 2 |l — A~ r |. The semiclassial contribution (42) 
comes from all the trajectory pairs for which one trajectory winds r more times around the periodic 
orbit p than its partner. The contribution comes from the close vicinity of the origin, i.e. from 
trajectory pairs for which the total number of windings around the periodic orbit is very large. The 
diagrammatic rules that we have encountered so far are the first three in table 1, and the total 
contribution to C(e) follows as 

C p > r (e) W ^ M 2 (f! ie) 2 -V cos {- l -rS p + \r^ v + f rT^j . (43) 
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Finally we perform the sum over the channels which gives a factor M(M + k—1) as in section 3. The 
result is 

C*v( £ ) « 2»e M M(M + 2 W ~ 1) ^,, r cos (-IrS, + + 0(e 2 ) . (44) 

In the case of systems without time-reversal symmetry (k = 1) this yields exactly the contribution of 
the r-th repetition of the trapped periodic orbit to the time delay. However, for systems with time- 
reversal symmetry the prefactor is slightly wrong. The reason for this is the same as why the diagonal 
approximation in section 3 did not give the correct mean time delay. There are other correlations 
between trajectories that have to be included. So far we have considered trajectory pairs, a and a', 
which differ only in the number of windings around the periodic orbit p. However, the orbit a can 
also have additional close self-encounters. As a consequence, a and a' can differ also in the way in 
which different links are connected in the encounter regions of the self-encounters. We will see in the 
next section that this leads to a variety of further correlations. 



4.2 Combinations of periodic orbit encounters and self-encounters 

Consider a trajectory pair with self-encounters. The number and types of the self-encounters are 
specified by the vector v. In addition, the trajectory pair can have an encounter with a periodic orbit. 
We consider first the situation where the periodic orbit encounter occurs during one of the (L + 1) 
links. The additional encounter divides a link into two parts and increases the total number of links 
to (L + 2). The joint probability density for this case is 

/' T 
dti . . . dtL+i — J? • (45) 

+ ^-^ncnLl^nc 

where the prime denotes again that all link times, including the remaining one ti+2, have to be 
positive (see discussion around equation (20)). With this probability density one finds again that the 
contribution to the correlation function C(e) factorizes and that it can be obtained by the diagrammatic 
rules that have been derived in the previous sections. They are summarized in the first three lines of 
table 1. After a summation over the channels one obtains 

C^{e) « (L + 1) M{ ^J~ l) N{v){-l) v 2ie»A p , r cos (-±rS p + f r^) + 0(e 2 ) . (46) 

There is a factor (L + 1), because the periodic orbit encounter can occur during any of the original 
(L + 1) links. 

One can also have the situation that a periodic orbit encounter overlaps with a self-encounter. 
In other words a self-encounter occurs in the vicinity of a periodic orbit. This leads to interesting 
consequences. In the simplest situation a two-encounter occurs in the vicinity of a periodic orbit in a 
system with time-reversal symmetry. A Poincare section for this case is shown in figure 4. An orbit 
a, indicated by the full lines, has a first encounter with the periodic orbit, and its intersection points 
with the Poincare section follow a hyperbola in the downward direction (the direction is indicated by 
arrows). Then the trajectory goes away from the periodic orbit, makes a loop and comes back and 
follows the periodic orbit in the opposite direction. The (time-reverses of) the intersection points with 
the Poincare surface then follow a second hyperbola in the upward direction. (When we speak of the 
intersection points of a in the following, we mean this to include the intersection points of its time 
reverse.) 

This situation is different from a usual two-encounter where a has only two intersection points 
with the Poincare surface, (ui,s±) and (ii2,s-2). In a usual two-encounter the partner orbit a' has 
intersection points which are approximately given by (u±,S2) and (u2,si), i.e. they are obtained by 
drawing a rectangle with the intersection points of a in opposite corners. If we apply this picture to find 
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Figure 4: An example of a trajectory (the dots on the two full lines) with a two-encounter in the vicinity 
of a periodic orbit and some of its partner trajectories (the dots on the dashed lines). The partner 
trajectories are obtained by drawing rectangular boxes which have points of the original trajectory 
(full lines) in opposite corners. 



the partner orbit in figure 4 we see that we have several possibilities to use this construction to find a 
partner orbit which traverses the loop in the opposite direction, because a has many intersection points. 
Some of these possibilities lead to the same partner trajectory. For example, in figure 4 the rectangles 
that are between the two full lines yield a partner trajectory that follows the two dashed hyperbolae 
that lie between the two full hyperbolae. If other intersection points of a are connected by rectangles, 
then other partner trajectories are obtained. Figure 4 shows one further example by the dashed 
hyperbolae that lie outside the full hyperbolae. But there are more possibilities that are not shown 
in the figure. The number k of rectangles that correspond to the same partner trajectory follow from 
the Poincare mapping (29) with |A„| = e XpTp . It has the form k = V In —\ — r + \ -m In —\ — r. 

*■ 1 o \ / I P I Apip maxj \Uj\ A p ip maxj \Sj\ 

This suggests the following definition of the encounter time 

i£ n CT c « kT p « -L In --^ — , (47) 

A p maxj | Si | x maxj \Uj\ 

which has a close similarity to (21). 

The occurrence of several partner orbits can be understood in the following way. A partner 
trajectory a' of a differs in the direction in which the loop between the two periodic orbit encounters 
is traversed. But it can also differ in the number of times it follows the periodic orbit before and 
after the loop, as long as the total number of traversals is the same as for a. For example if a has 
k\ intersections in the square of the Poincare section with side lengths c before the loop, and &2 
intersections after the loop, then the partner trajectory can have k\ + d intersections before the loop 
and k2 — d intersections after the loop, where d is any integer with k\ + d > and — d > 0. This 
means that the presence of the periodic orbit leads to a large number of possible partner trajectories. 

However, we are interested in obtaining the periodic orbit contributions to the time delay. Hence we 
have to consider trajectories a" which have altogether r more periodic traversals than the trajectories 
a' . These trajectories a" then differ from the original trajectories a in the number of periodic orbit 
traversals as well as in the direction in which the loop is traversed. In the following we calculate the 
semiclassical contribution of the trajectory pairs a and a". 

Let us consider the general case that an /-encounter occurs near a periodic orbit p. A partner orbit 
a' is obtained by reconnecting the links in a different way in the encounter region. This reconnection 
is specified by a permutation n(j), j = 1, ... ,1 which consists of one cycle of length I [6, 25]. The 
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contribution of each link M(i—ie) 
contribution of each ^-encounter — M(l — He) 

r-th contribution of a periodic-orbit encounter 2ie^A Pir cos (— ^rSp + f + ^e/irT p ) 

r-th contribution of periodic-orbit plus /-encounter 2\le\xA p ^ r cos (—-^rS p + ^rfi p + ^efj,rT p ) 

Table 1: Diagrammatic rules for the different contributions to the correlation function C(e). 

action difference in [6, 25] can be expressed in the form AS* = X^=i s j( u j ~ u ^(j))- ^ we now consider 
partner trajectories a" that have additional periodic orbit traversals r±, . . . , r\ during the I encounters 
with the periodic orbit, with Y2j r j = r > then the action difference with a changes according to (31) 
and is given by 

l 

AS = Sj(uj - u n{j) Kp r ' 3 )-rS p = s T Bu - rS p (48) 

3=1 

where s and u are the vectors of the I coordinates in the stable and unstable directions, and B = I + B 
where Bji = —8^fj\h. p rj . We have det B = 1 + det B, because the permutation consists of one cycle of 
length I. Furthermore, (I — 1) column exchanges bring B into diagonal form, and hence det B = — A£, 
and det B = 1 - A p . 

If we have one l-encounter a in the vicinity of a periodic orbit p, and no further self-encounters, 
then the probability density is given by 

/' T 
dti . . . dt t . (49) 

where the encounter time is given in (47). The factor 1/k = T p /tenc takes care of the overcounting of 
trajectory pairs. The contribution to C(e) follows from (36), (37) and (49) as 

E ^ ~ T i r - / ds du wj* II / d ^ ^ {1 - le)u + - ( 5 °) 

The expression factorizes again, and the contribution from the encounter region is obtained by ex- 
panding the exponential which contains the encounter time. The contribution originates again from 
the linear term of this expansion. The resulting diagrammatic rule is listed in the forth line of table 
1. Note that the contribution does not depend on how the repetition number r is split into parts 
ri, . . . , 77. In fact, the integral over the s and u coordinates sums over the different ways of splitting r 
into parts, because a trajectory a has many intersection points. The different s and u coordinates for 
these intersection points correspond to different ways of splitting r, as in figure 4 and the discussion 
after equation (47). 

Finally, we include the possibility that there are additional self-encounters which are not near a 
periodic orbit. One finds that also in this case the contributions factorize, and they can be evaluated 
by the rules in table 1. If the total number of self-encounters is given by the vector v, and we sum 
over all the cases where one of the V self-encounters is in the vicinity of the periodic orbit p, we obtain 



C^ r (e) « f-I>) M( mL + J + 2 l) N{v){-l) v ^A p , C os (-±rS p + |r^) + 0(e 2 ) . 



(51) 
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This expression differs from (46) in that the factor (L + 1) is replaced by — J2 a la = When we add 
the two this factor becomes one. Finally we sum over all vectors v and include also the contributions 
from section 4.1 (corresponding to L = V = and N(v) = 1). 

oo L—V+l=k pj, \/_-i\V / 1 \ 

C^(e) « (M + k - 1) £ E Mk 2ie » A P>r cos ( -jfS, + |r Mp + 0(e 2 ) 

fc=l u ^ ' 



2ie/j,A Pjr cos ^-irS'p + ^ r ^pJ + 0(e 2 ) . 



(52) 



where the sum rule (27) for N(v) has been used. 

Equation (52) is the final result of this section. Using (10), it gives the correct contribution of the 
r-th repetition of the periodic orbit p to the time delay (3) in systems with or without time reversal 
symmetry. 

We have seen in this section that the vicinity of a periodic orbit gives rise to a rich variety of 
correlations between trajectories. We concentrated on those correlations that are responsible for the 
semiclassical periodic orbit terms of the time delay. There are further correlations which we do not 
explore in the present article. One example are the trajectory pairs of this section with r = 0. There 
can also be multiple encounters with a periodic orbit p which do not correspond to the picture that 
one of the usual self-encounters occurs near a periodic orbit. An example is a trajectory which visits a 
periodic orbit twice and follows it both times in the same direction and has no further self-encounters. 
Furthermore, a trajectory can have encounters with several different periodic orbits. These terms, 
however, appear in higher order terms in the expansion in e. (Visits of n periodic orbits appear in 
terms of order e n .) These additional correlations merit further study. 

An interesting question is whether one obtains periodic orbit contributions to the conductance 
with the type of open orbit correlations in this section. The Landauer-Biittiker formulation for the 
conductance can be obtained, in certain situations, from linear response theory (Kubo formula). On 
the other hand, when the Kubo approach is applied to the bulk conductivity in antidot lattices it gives 
semiclassical expressions in terms of periodic orbits [24]. It is an open question whether periodic orbit 
contributions exist for the Landauer-Biittiker conductance as well. 

The conductance is proportional to the total transmission through the cavity 

T= E S ^ ai AE)S* bouta jE) , (53) 

where the only difference to C(e) lies in the different channel sum (the sum is over incoming and 
outgoing channels), and e = 0. Since the periodic orbit contributions to C(e) in (52) vanish in the 
limit e = one finds that there is no periodic orbit contribution to the conductance as well in this 
leading order semiclassical calculation. This does not completely rule out the possibility of periodic 
orbit terms in the conductance, but if they exist they need to have a different sign in the reflectance, 
because transmission and reflectance have to add up to a constant, the number of channels in the 
incoming lead(s). 



5 Correlation functions of the time delay 

In this section we look at a correlation function of the time delay, and show that it is possible to get 
the leading order result of RMT using open trajectories. This addresses the concern of [15] that the 
diagonal approximation in the approach using trajectories gave a different result. The key is that there 
are other contributions of the same order as the diagonal approximation for open trajectories, and 
these need to be included to get the leading order term. The calculation for the correlation function 
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of the time delay is similar to that for the conductance variance and, in particular, to the Ericson 
fluctuations [11]. For this reason we shall be brief and only sketch the calculations. The correlation 
function that we consider is defined as 



R 2 (u;,M) = ^ A " d)/E " . (54) 



Again it is convenient to specify the energy difference in units of M(2ird)~ 1 . The correlation function 
involves an energy average over an energy range AE which is classically small, but encompasses many 
resonances E » AE >> l/d(E). In RMT the leading order result for a large number of channels 
M is given for the two considered symmetry classes by [12, 13] 

When we express the correlation function in (54) as a sum over quadruples of trajectories we subtract 
the term f^, by removing the trajectories pairs of section 3 that give the mean delay time. Inserting 
the semiclassical approximation for the time delay (6) results in 



R 2 ( UJ ,M) = ^(Y j T a r^A a ^*,^^,e^ 5 «- s «' + ^-V)e^( T «- T «'- r /5+ T /3')\ . (56) 

H \ a,b a,a'{a->b) I 

c,d p,fl>(c^d) E 

Here, a, a' are trajectories from channel a to b and (3, [3 1 are trajectories from channel c to d and we 
sum over trajectories and channels. The prime at the second sum indicates that we have removed the 
orbits where a a' and (3 « (3'. 

We start with the calculation of the diagonal contribution from [15]. Because we have removed the 
trajectory pairs that give the average time delay, we have to consider only the cases when S a = Sp> 
and S/3 = S a r. Without time reversal symmetry this means that a = (3' and (3 = a', which requires 
that the channels a = c and b = d. Hence we get a factor M 2 from the sum over the channels. With 
time reversal symmetry there are three additional cases. We can also have a = f3' and (3 = a', where 
the over line indicates the time reverse. This requires a = d and b = c and the channel sum gives M 2 . 
Furthermore, we can have a = (3' and (3 = a', or a = (3 1 and (3 = a' . Both cases require a = b = c = d 
and each channel sum yields M. The total channel factor for the systems with time reversal symmetry 
is thus 2M(M + 1), and in leading order for large M the channel factor for the two symmetry cases 
is ~ kM 2 . We can perform now the sums over the trajectory pairs by replacing them by an integral 
according to (13), and we perform the sum over the channels by multiplying with the channel factor 

£f*(c;,M)«^- / / dT^r^e-^+^e 1 ^-^), (57) 

Jo Jo 

which, after integrating, is 



We have a factor of two different from the expected result (55), and a different functional form. To 
correct for this we have to add sums over quadruplets with encounters that contribute at the same 
order of M _1 . The reason why they can contribute at the same order as the diagonal term is because 
they have larger channel factors. 

One example is a quadruplet where the trajectories a and (3 have one close encounter. Then their 
partner trajectories follow one of these two trajectories from the entrance channel to the encounter 
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Figure 5: An example of a trajectory quadruplet where the trajectories a and j3 (full lines) have one 
two-encounter. Their partner trajectories (dashed lines) cross the encounter region differently. The 
encounter region is indicated by a rectangular box. 



region where they switch over to the other trajectory and follow it to the exit channel, see figure 5 
(or figure 4j in [11]). If we label the link time betweens entrance or exit channels and the encounter 
region by ij, then the times along a and (3 can be written as 

T a = t\+ t enc + t 2 , Tp = i 3 + t enc + t 4 . (59) 

If a' follows first a and then /3, and (5' does vice versa, we have 

T a ' =t\ + tenc + *4 , Tpi = t 3 + t enc + i 2 • (60) 

These quadruplets are only possible if channels b and d are identical. We can also interchange the role 
of a' and (3' which requires that a = c. In systems with time reversal symmetry we have to consider 
also cases where the role of a' and/or is taken over by its time reverse, similarly as for the diagonal 
approximation. The sum over trajectories for these different configurations always gives the same 
result, hence it is sufficient to do the calculation for the case (60) and to perform the channel sum by 
multiplying with the total channel factor. In systems without time reversal symmetry this channel 
factor is 2Af 3 , and in system with time reversal symmetry it is 4M 3 + 8M 2 + AM. In leading order it 
is hence ~ 2kM 3 . 

The difference in times for the case (60) is 

-{Ta-Tof -Tp-Tp') = t 2 -t 4 , (61) 
and the contribution to the correlation function is given by 

flte 



rp4 



f dsdu (ti + i enc + t 2 )(t 3 + i enc + U)^- — et su . (62) 

J ^ '*enc 



This then can be evaluated with the rule (23) and yields 

2k uj 2 
M 2 (1 + w 2 ) 2 



(63) 



Following a similar process we find the contribution from the other diagrams in figure 4 of [11] and give 
the results in table 2. When we now multiply these contributions by their channel factor M 4 + M 2 ~ 
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4c 

-l-2o; 2 


4d 

8u; 2 (l+u 2 ) 


4e 

2-2w 2 -8a; 2 (lW 2 ) 


M b (l+w 2 ) 2 


M b (l+uj 2 ) 2 


M b (l+w 2 ) 2 



Table 2: Contribution of different types of trajectory quadruplets to the time delay correlation function. 
The labelling corresponds to figure 4 in [11]. 



M 4 for the unitary case and add the terms calculated above, we get the following result of 

2 l-io 2 



M 2 (1 + w 2 ) 2 ' 



(64) 



which is indeed the leading order term of equation (55). For the orthogonal case, the channel factor 
for the configurations in table 2 is M 4 + 2M 3 + 3M 2 + 2M ~ M 4 to leading order. However there are 
additional trajectory quadruplets that contribute (4f-4h) which are related to the quadruplets 4c-4e 
by time reversal of parts of the structure. In fact these additional quadruplets just give the same 
contribution again as those in figures 4c-4e. Therefore, the leading order contribution is simply twice 
that for the unitary case, again in line with equation (55). 

It is worth noting here that using the symmetrized version of the time delay, equation (10) with 
(9), we get a different result for each trajectory quadruplet, but the sum of their contributions gives 
the same result as here. To calculate higher order terms, using the methods of [11], the similarities 
can be exploited by defining a 'symmetrized' correlation function which can be written in terms of the 
correlation function of the scattering matrix elements. 



-1 d 2 
M 2 deide 2 



C a [e 1 ,E + 



ujM 
And 



£1=62=0 



(65) 



where C fl (e, E) is the semiclassical approximation to the fluctuating part of the correlation function 
of the scattering matrix elements at energy E. This means that the trajectory pairs responsible for 
the average part are removed. The calculation is very similar to that of the Ericson fluctuations. A 
complication in comparison to the calculation of the mean time delay in chapter 3 is that one does not 
have one simple diagrammatic rule for every link. Each link is traversed by two of the four trajectories 
a, a', P, (3' , and the diagrammatic rule depends on which two of the four trajectories that are. The 
summation over the contributions from the different families is best done by a computer. 
A final point is that if we consider C fl (0,£) C R (0,E) we are calculating 



(Tr S(E)S^E) -m) 



(66) 



which should be zero because of the unitarity of the scattering matrix. We checked this by using the 
formulae for the conductance variance [11], but with the appropriate channel factors that are given 
in this section. The result is indeed 0, and the unitarity of the scattering matrix is preserved by the 
semiclassical approximation if all semiclassical contributions are included. 



6 Conclusions 

The time delay is an interesting quantity to study because of the two alternative semiclassical de- 
scriptions of it. The picture in terms of scattering trajectories is similar to the conductance, and by 
using the semiclassical methods of [11] which includes trajectories with self-encounters we obtained 
the average time delay. We considered also a correlation function of the time delay. Here the diago- 
nal approximation is not enough to obtain the leading order term (as noted by [15]), because more 
complicated trajectory quadruplets with encounters contribute at the same order. By including these 
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contributions we find that the unitarity of the semiclassical scattering matrix is restored, and that the 
semiclassical method does indeed provide an accurate description in agreement with RMT. This is as 
expected from other correlation functions in [11]. 

The main result of this article, however, is the discovery of a new type of scattering trajectory 
correlations that recreate the periodic orbit terms of the time delay. These trajectories approach a 
periodic orbit, follow it closely some number of times, and then part from it. Approaching the orbit 
guarantees that they remain in the system, as long as the periodic orbit does and they remain close to 
it. This is the reason why only the trapped periodic orbits appear in the orbit sum of the time delay. 
For systems without time-reversal symmetry these periodic orbit encounters are enough to recreate 
the periodic orbit terms of the time delay, but for systems with time-reversal symmetry a small 
additional contribution is needed. We obtained this contribution from combinations of self-encounters 
and periodic orbit encounters. 

We found that the occurrence of self-encounters in the close vicinity of a periodic orbit leads to 
a wealth of new possible correlations between trajectories. We considered in this article only those 
correlations that are needed for the periodic orbit terms of the time delay, but there are further possible 
correlations. These new type of correlations deserve further study, and it can be expected that they 
play a role also in other contexts. For example, they might be relevant for periodic orbit correlations 
in closed systems, as envisioned in [29]. For the particular case of the Landauer-Biittiker conductance, 
however, we found that the trajectory correlations that yield the periodic orbit terms for the time 
delay do not give similar periodic orbit terms for the conductance. 
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